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SUMMARY 


An alternative method of analysis to determine the dynamic response of 
structures subjected to base accelerations is presented. The method is exact as 
opposed to the approximate technique of using unusually large masses and loads to 
enforce desired base accelerations. This paper presents the relevant equations of 
motion, ALTERS for direct and modal frequency-, random- and transient-response rigid 
formats, and illustrative examples. 


INTRODUCTION 


Dynamic environments of industrial structures and structural components are 
frequently specified in terms of base accelerations. The dynamics analysis capa- 
bilities of NASTRAN, however, provide for the specification of loads rather than 
accelerations. An approximate technique is therefore, generally practiced wherein 
the structural degrees of freedom with known (base) accelerations are assigned very 
large masses (or mass moments of inertias, of the order of 10^ to 10^^ times the 
total structural mass) and subjected to correspondingly large loads to enforce the 
desired base accelerations (ref. 1). This method can, in some instances, lead to 
erroneous results as shown by the frequency response function in Figure 1. The 
function represents 0^^^^ in element 2 when the cantilevered plate, shown in 
Figure 2, is subjected to a unit base acceleration z (f) =1. A concentrated mass 
of 10^ units was used in the z (translation) degree of freedom at grid point 1. The 
solution was obtained using the direct frequency and random-response rigid format 
DISP RF 8. A viscous damping matrix proportional to the structural stiffness matrix, 
BDD = 2.0 E-6 KDD, was used enabling a direct comparison of the results with those 
from the alternative method. 

The method discussed in the following sections avoids the use of fictitious 
large masses thereby eliminating any associated conditioning problems. For a modal 
formulation of the problem, the modal basis currently available in NASTRAN rigid 
formats is used. These modes are the base-relaxed modes including rigid body modes, 
and constitute a plausible basis for base excited dynamic response calculations. 

The prevailing boundary conditions are satisfied by all modes. The user directly 
specifies base accelerations on existing NASTRAN bulk data cards. 


NOMENCLATURE 


Note: A consistent set of FLT units has been used throughout this paper. 
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Damping matrix, FL 
Viscous damping, FL 
Frequency, T ^ 

Modal damping coefficient 

Stiffness matrix FL ^ 
Stiffness, FL ^ 

-1 2 

Mass matrix, FL T 
-1 2 

Mass, FL T 

Load vector, F 

Derivative operator ~ 

dt 
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Time, T 

Displacement vector, L 
Translational accelerations, LT 
Circular frequency = 2Trf, T 
Circular frequency of 1th natural mode, T 
Modal matrix 

Normal stress component in basic x direction, z- fiber, FL 
Subscripts: 
a, d, h, Z, p, r 
D = d-r set 


-2 


NASTRAN displacement sets 


METHOD OF ANALYSIS 


Direct Formulation 

After the application of constraints and partitioning to both structural and 
direct input matrices, the equation of forced motion is 
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The displacement vector Is partitioned as 


{u^} “ {Pj (t or ^)} . 
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with Up representing the base accelerated degrees of freedom. Equation (1) can now 
be rearranged Into the following two equations: 
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2 2 
Given the base accelerations p Up, equation (3) is solved for u^, pu^ and p u^. 

Equation (4), in turn, can be solved for Pp — the loads required on the base degrees 

of freedom to cause the desired base accelerations. 


Modal Formulation 

The displacement vector uj is written as 


{Ud> = 


II 

■‘^Dh 

lu,/ 

-^rh- 




(5) 


where the modal matrix has been appropriately expanded to include any extra 
points (ref. 2). Substituting equation (5) in equation (3), and premultiplying both 
sides by the resulting equation of motion is 


where 
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(continued) 




91 





T 

I 


0 





■ i:*dhi 


[M ]p^{u } 


ORIGINAL PAGE 19 

OF POOR QUALITY (continued) 


(7) 


Consistent with the existing capabilities of NASTRAN rigid fornats for modal 
frequency- and random-, and transient-response analyses (DISP RFs 11 and 12), the 
modal damping matrix consists of contributions due to, 

1. the damping matcix proportional to the mass matrix as 



I 


with the elements of the diagonal matrix corresponding to all extra points set to 
zero, and 


2. the direct input matrix. 

2 2 
Gi^'en the base accelerations p u^., equation (6) is solved for uj^, puj^ and p u^. 

Equation (5) is used to recover the displacement vector uj and its rates pu^j and 

p^uj. Equation (4) can be used to recover the loads Py The damping terms B^p and 

are partitions of the directly specified damping matrix B2PP, 


IMPLEMEOTATIOK IN NASTRAN 


The method of analysis discussed in the previous section has been implemented 
in NASTRAN April 1982 release in the form of DMAP ALTERS. The ALTER packages for the 
displacement approach rigid formats 8, 9, 11 and 12 are given in the Appendix. 
Existing NASTRAN utility modules have been used to partition and merge various 
matrices for the rearranged equation of motion. The functional module FRRD2 for the 
frequency response rigid formats has been modified to solve coupled equations of 
motion. These modifications are also included in the Appendix. 


In using these ALTER packages, the following points are to be considered. 

1. The base accelerated degrees of freedom are specified on the SUPORT bulk data 
card. 

2. The base accelerations are specified on RLOADi or TLOADi bulk data cards akin 
to specifying loads. 

3. The base accelerated degrees of freedom must have non-zero mass (or mass moment 
of inertia). No fictitious large masses are required. 

4. In rigid formats 8 and 11, external loads can be applied to all (p-r) set 
degrees of freedom. 

5. In rigid formats 9 and 12, external loads can be applied to all (d-r) set 
degrees of freedom. 

6. An OLOAD request for the base accelerated (r set) degrees of freedom in RFs 8 
and 11 results in the Joads on these degrees of freedom necessary to cause the 
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specified base accelerations. Such a request in RFs 9 and 12 will output the 

specified base accelerations. This is due to the non-availability of puj and 

p^Ujj in distinct data blocks. 

7. If OLOADs are requested for the base accelerated degrees of freedom in RFS and 

11, a real diagonal matrix named FREQ must be input on DMI bulk data cards with 
entries sorted in an ascending order from the FREQ or FREQi bulk data card. 

The order of the FREQ matrix must be equal to the entries on the FREQ or FREQi 

bulk data card. 

8. Mode acceleration method of data recovery is available both in RFll and 
RF12. 

9. The data recovery procedures in all the four rigid formats remain unchanged, 
with the exception of stress recovery. The stresses are computed using dis- 
placements (or modal displacements) relative to the base. This is due to the 
fact that in problems with specified base accelerations , the absolute displace - 
ments can become extremely large. The subsequent stress calculations, as a 
result, are based on small differences of large numbers, and can be in error as 
shown by the stress response at very low frequencies in Figure 1. It is to be 
noted that any limitations imposed by these ALTER packages are as a direct 
result of utilizing existing functional and utility modules, with the necessary 
exception of FRRD2. These limitations can easily be overcome by creating new 
(dummy) modules. 


ILLUSTRATIVE EXAMPLES 


Two problems (Figures 2 and 3) are considered to illustrate the alternative 
method of determliing base accelerated dynamic response. The problem in Figure 3 is 
used to illustrate the accuracy of all the four ALTER packages. The problem in 
Figure 2 is used to compare response calculations with those shown in Figure 1. 


Figure 3 shows a 2-degree of freedom system subjected to a known base accelera- 
tion at mass m^^. The problem is to determine the acceleration response of mass m 9 . 
The following steps are followed in obtaining and cross-checking the solution by 
various ways: 


1. RFS with ALTERS is used to determine y 2 (f) and Pj^(f), given yj^(f). 

2. RFS without ALTERS is used to determine yi(f) and y 2 (f), given Pj^(f). 

3. Steps 1 and 2 are repeated with RFll. 

4. RF9 with ALTERs is used to determine y 2 (t), given yj^(t). 

5. Step 4 is repeated with RF12. 

Figures 4 and 5 show the results of steps 1 through 3 which compare well with 
the theoretical results given by 


y 2(‘^0 


k + Icoc 


k + io)C - o)^m„ 


y^(w) ,and Pj^(tJ) = mj^yj^(i-j) +m^y^{.) 


(9) 
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y„(t) = 1 : cos (v^l - zj a) t + Y.') , 

2 /r~ r2 n o 


where 


4 = c / (2/km ) , 


u> = /k/m^ , and 

n 2 


Yq" = tan ^ ( ^ ) , 


/I - 4 2 


( 10 ) 


for 


y^(t) = 0 , t < 0 


= 1 , t > 0 . 


Figure 7 presents the frequency response function of Figure 1, using RF8 with 
base acceleration ALTERS. For comparison, a damping matrix proportional to the 
structural stiffness matrix was used in both solutions (BDD =2.0 E-6 KDD) . The 
response does not rely on the selection of any fictitious large masses. The error in 
the stress response at low frequencies is also eliminated. 


CONCLUDING REMARKS 


1. An alternative method of determining base accelerated dynamic response in 
NASTRAN has been presented and demonstrated, avoiding the use of fictitious large 
masses . 


2. Although this paper discusses the problem and its solution in terms of absolute 
degrees of freedom, a number of variations can be simply achieved to suit particular 
problems. As an example, the introduction of the degrees of freedom relative to the 
base degrees of freedom, at least in the absence of extra points, can be easily 
accomplished by the transformation 


with 



I ! Dl/^£j^rel.\ 

5TiJ\"u™7 


[D] - [KjJ . 


(set e null) 




( 11 ) 


This is useful in addressing shock spectrura response problems, and can lead to 
symmetric coefficient matrices and uncoupled (modal) equations of motion. 


4. The method can form the basis fov considering displacement and velocity base 
excitation problems in NASTRAN. 
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APPENDIX 


$$ 

$$ ALTERS TO SOL 8 (NASTRAN APRIL 1982) FOR BASF AIXFLERATI ON . 

«$ 

ALTER 87 $ 

RBMGl USET»KAAr / KLL f KLR » KRR » f » * 

RBMG2 KLL / LLL $ 

RBMG3 LLL.fKLR»KRR / MM $ 

ALTER 12.1. >122 $ 

FRLG LASEXX > USETM > MLT r F RL » GMD » GOM > MI T » /PPF » PSF PDF » FOl. » PHF/ 

♦DIRECr*/FREOY/*F REO» % 

EQUIO PPF»PMF/NOBET $ 

9EC USETM/OP/»P*/)K(:;OMP#/)|tR)»t $ 

OEC USETM/OM/*M*/)KC;OMP)K/#R* $ 

PAR I N PPF » > OP/PPF 1 » PPF 2 > > /I * 

PARTN PMF » > OM/ PDF 1 » PMF2 > > / 1 * 

PARTN MMM > VM > /MMM 1 1 > MMM2 1 . MMM 1 2 v MMM22 $ 

PARTN BMMfOMr/BMMll »BMM21fBi:iM12>BMM22 % 

PARTN KMM > OM > /KMMl 1 > KDM21 > KMMl 2 > KMM22 $ 

MERGE MMDll > » >Mi:iM22>yMv/ MMMl /-l/O/l ♦ 

MERGE BMM 1 1 > > BMM 1 2 » » VM » / BMM 1 / - 1 70/ 1 $ 

MERGE KMMll » >KMM12> >DM>/ KMMl /-l/O/l $ 

MPYAM MMM12>PPF2>PMF1/ PMFF 1 /0/-1 ♦ 

MPYAM MMM22»PPF2» / PMFF2 /O $ 

MERGE PMFF1»PMFF2> > » > 01 . 1 / PMFF /1/0/2 $ 

F RRM2 KMMl > BMMl > hMMl > > PMFF > F OL/UMOF /0 .0/0 . 0/-1 . 0 $ 

AMM FRE0>/0MEGAI/(0.0>6.28318S) $ 

ADD FRE0>/0MEGA1/(6.283185>0.0) $ 

DIAGONAL 0MEGAl/0MEGA2/*SQUARE*/2.0 ♦ 

MERGE MMM21 > >MMM22> >UM> / MMM2 /1/0/2 $ 

MERGE BMM21 > >BMM22> >0M> / BMM2 /I/O/ 2 ♦ 

MERGE KMD21> >KMM22> >0D> / KMM2 / 1/0/2 ♦ 

SMPYAD MMM2>UM0F>0MEGA2> > > / TEMPI / 3/-1 $ 

SMPYAM BMM2>UM0F>0MEGAI>>> / TEMP2 / 3 ♦ 

ADD TEMPI >TEMP2 / TEMP3 ♦ 

MPYAM KMM2>I.JM0F>TEMP3 / PPF2/ 0 % 

MERGE PPF1>PPF2> > > >0P / PPF / 1/0/2 ♦ 

EQUIO PPF>PNEF/MPCF1 $ 

CONM LBL18A>MPCF1 $ 

VEC USETM/OECNEM/*P#/»NE*/»M* $ 

PARTN PPF > >OFCNEM/PBNEF>PMF>> / 1 $ 

MPYAM GMM>PMF>PBNEF/ PNEF / 1 $ 

LABEL LBL18A ♦ 

EQiJio pnef>pfe:f/single. ♦ 

CONM LBi..l8B> SINGLE * 

OEC USETM/OE;CFE:S/*NE.)|£/#FE*/)KSi|t $ 
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PARTN PNEF»»VKCFES/PFEF»PSF.» QUALITY 

LABEL LBL18B $ 

EQUIV PFEF» PDF/OMIT $ 

COND LBLlSCrOMlT * 

VEC lJSETi:i/VECBO/*FE*/)|tD*/)KO* ♦ 

PARTN PFEF»rVECi;it)/PBDFfPOF»» / 1 ♦ 

MPYAD BOD»POF»PBDF/ PDF/ 1 ♦ 

LABEL LBL18C $ 

ALTER 123^123 $ USE FOL INSTEAD OF PPF 

VDR CASEXX » EOD YN » USETD » UDVF » FOL » X YCDB r /OUDVCl r /*FREQRESP*/ 

*DIRECT*/S » N » N0S0RT2/S » N » NOD/S » N » h'OP/0 * 

ALTER i39fl39 $ USE FOL INSTEAD OF PPF 
SDR2 


UEC 

UEC 

PARTN 

PARTN 

MPYAD 

UMERGE 

UMERGE 

EOUIU 

COND 

SDRl 

LABEL 

SDR2 


CASEXX » CSTM» MPT »DIT*EODYN»SILDf » »BGPDP» FOL »QPC»UPVC» EST»X YCDB » 
PPF/OPPCl fOQPCl fOUPUCl f »OEFCl »PUPUCl/*FREORESP*/ 

S»N>N0S0RT2 « 

USETD / UD2 / $D$/*A*/*E$ $ 

USETD / UA2 / $A$/$L$/*R* $ 

UDUF»»UD2 / UAUFfUEVFff / 1 ♦ 

UAUF»»UA2 / ULVF»URVF»» / 1 $ 

DM»URVFrULVF / ULVFP / 0 / -1 * 

USETD fULUFPf / UAUFP / #A*/*L*/*R* ♦ 

USETD »UAVFP»UEVF / UDUFP / *D*/*A*/#E.1t $ 

udufp»upvcp/noa « 

LBL19ArN0A t 

USETD»»UDUFP»»»GOD»GMD»»»» / UPUCPrr / 1/*DYNAMICS* $ 

LBL19A $ 

CASEXX f CSTM » MPT » DI T » EOD YN r SILD r r » BGPDP » FOL » r UPVCP » ES T » 
XYCDB»/»»»OESCl»» /»FREQRESP»/SfNrN0S0RT2 ♦ 


ENDALTER 
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♦t ALTERS TO SOL 9 (NASTRAN APRIL 1982) FOR BASE ACCELERATION, 

«$ 

ALTER 87 t 

RBM61 USETfKAA* / KLL»KLR»KRRf » » » 

RBHG2 KLL / LLL $ 

RBHG3 LLL»KLR>KRR / DM $ 

ALTER 125fl25 * 

9EC USETD/VD/*D*/*COMP*/#RiH i 

PARTN PD»»Vn/PDl»PD2f »/l $ 

PARTN HD0f0D>/MDDllfMDD21>MDD12»MDD22 $ 

PARTN BDD»VDf/BDDll»BDD21»BDD12rBDD22 $ 

PARTN KDD»9D»/KDD11»KDD21fKDD12tKDD22 $ 

MERGE MDD11» » »MDD22»VD»/ MDDl /-l/O/l ♦ 

MERGE BDDll»»BDD12»fVD»/ BDDl /-l/O/l » 

MERGE KDDll»fKDD12, .VD»/ KDDI /-l/O/l » 

MPYAD MDD12»PD2»PD1/ PDDl /0/-1 $ 

MPYAD MDD22»PD2» / PDD2 /O ♦ 

MERGE PDDl»PDD2ip»»fOD/ PDD / 1/0/2 ♦ 

TRD CASEXX » TRL f NLFT f DI T » KDDI » BDDl » MDDl » PDD/ 

UDOTrPNLD / *DIRECT*/NOUE/l/S»Nf NCOL/C» Y» ISTART $ 

ALTER 139»139 ♦ 

CASEXXfCSTM»MPT»DIT»EQDYN»SILDf » »BGPDP»TOL»QPf UPO»ESTf 
XYCDB»PPT/0PP1 »0QP1 fOUPVl r rOEFl »PUGO/*TRANRESP* $ 

USETD / VD2 / *D»/»A*/*E» ♦ 

USETD / VA2 / #A*/»L#/*R» ♦ 

UD0T»fVD2 / UAVT,UEVT»f / 1 ♦ 

UA0T»»VA2 / ULVTfURVTff / 1 ♦ 

DMfUROTfULVT / ULVTP / 0 / -1 ♦ 

USETD »ULVTP» / UAOTP / *AHt/*L#/*R# ♦ 

USETD rUAVTP»UEOT / UDVTP / !KD#/*A#/#E)|c $ 
UDOTP»UPVTP/NOA $ 

LBL19ArN0A ♦ 

USETDf »UDOTP»f »GOD»GMD» / UPVTPrr / 1/*DYNAMICS* $ 
LBL19A $ 

CASEXXfCSTM>MPT»DITfEQDYN»SILD» r ,BGF'IiP»TOLf »UPVTP»EST» 
XYCDB»/f »»0ES1»» /*TRANRESP* $ 

ENDALTER 


SDR2 

VEC 

VEC 

PARTN 

PARTN 

MPYAD 

UMERGE 

UMERGE 

EQUIO 

COND 

SDRl 

LABEL 

SDR2 
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♦♦ ALTERS TO SOL 11 (NASTRAN APRIL 1982) FOR BASE ACCELERATION* 

«« 

ALTER 93 t 

PARAM //«A0D»/KDEK2/N0GENL/N0SIMP $ 

ALTER 95f96 « 

PARAH //*ADD*/NOBGG/-l/0 ♦ 

PARAH //*ADD»/N0K4GG/-l/0 $ 

GKAD USETD » GM » GO f K AA » » MAA f » K2PP » H2PP » B2PP/ 

KDD f BDD f NDD » GMD f GOB t K2D0 1 M20D » B20D/»FREORESP«/»D1 SP«/ 
*DIRECT*/0 . 0/0 ♦ 0/0 . 0/N0K2PP/N0H2PP/N0B2PP/HPCF1 /S INGLE/ 
ONI T/N0UE/N0K4GG/N0BGG/KDEK2/0 » Y » MODACC-- 1 ♦ 

EQUIV B2DD»BDD/N0BGG/H2DD»HDD/N0SIMP/K2Dn»KDD/KDEK2 $ 

FRLG CASEXX » USETD» DLT f FRL f GMD f GOD f DI T f /PPF f PSF f PDF f FOL f PHF/ 

»DIRECT*/FREQY/»FREQ» ♦ 

EQUIV PPFf PDF/NOSET ♦ 

VEC USETD/VP/*P*/*COHP*/*R» ♦ 

VEC USETD/VD/»D*/»COHP)|c/»R» $ 

PARTN PPF? fVP/PPF1fPPF2f f/1 ♦ 

PARfN PDFf fVD/PDF1fPDF2f f/1 ♦ 

PARTN MfjriFVDF/HDDllFMDD2lFMDD12FMDD22 $ 

PARTN BDDfVDf/BDD11fBDD21fBDD12fBDD22 $ 

PARTN KDDfVDf/KDD11fKDD21fKDD12fKDD22 ♦ 

MERGE MDD11fffMDD22fVDf/ MDDI /- VO/1 * 

MERGE BDD11ffBDD12ffVDf/ BDDl /-l/O/l $ 

MERGE KDD11ffKDD12ffVDf/ KDDI /-l/O/l $ 

MPYAD MDD12fPPF2fPDF 1/ PDFFl /0/-1 $ 

MPYAD MDD22fPPF2f / PDFF2 /O $ 

MERGE PDFFl F PDFF2 f F F fVD/ PDFF /1/0/2 ♦ 

GKAM USETDfPHIAfMIfLAMAfDITf f f f CASEXX/ 

MXHH F BXHH f KXHH f PHI DH/NOUE/C f Y f LM0DES=0/ 
CfYfLFREQ=0.0/CfYfHFREQ=-1,0/-1/-1/-1/ 
SfNfNONCUP/SfNfFMODE ♦ 

MPYAD PHIDHfPDFFf/PHF/1/1/0/0 $ 

SMPYAD PHIDHfMDDI fPHIDHf F f/MHH/3/I/j./O/I » 

SMPYAD PHIDHfKDDIfPHIDHfff/KHH/S/I/I/O/I $ 

PURGE B2HH/N0B2PP $ 

COND LBL13AfN0B2PP $ 

SMPYAD PHIDHfBDD1fPHIDHfff/B2HH/3/1/1/0/1 $ 

LABEL LBL13A $ 

EQUIV MAAfMIDD/NOUE * 

COND LBL13BfN0UE $ 

VEC USETD/VD1/*D*/»A*/*E* ♦ 

MERGE MAAffffVDIf/ MIUD /-l/O/l i 

LABEL LBL13B $ 

PARTN MlUDFVDF/MlDDllFrMlDD12F $ 

MERGE M1DD11ffM1DD12ffVDf/ M1DD2 /-l/O/l $ 

SMPYAD PHIDHfM1DD2fPHIDHfff/M1HH2/3/1/1/0/1 ♦ 

DIAGONAL MXHH/MINV/)»tSQUARE*/-l .0 ♦ 
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SMPYAD MINV»BXHH»M1HH2» » » / TEMPI / 3 $ 

ADD TEMPI »B2HH / BHH $ 

PARAM //*ADD*/NONCUP/l/0 ♦ 

ALTER 99 » 100 $ 

FRRD2 KHHfBHH»MHHt,PHFfF0L/UHVF/0*0/0.0/-1.0 $ 

DDRl UHVF»PHIDH/UDVF $ 

ADD FREQf/0MEGAl/(0.0f6«283185) « 

ADD FREQ»/0MEGAl/(6. 283185*0.0) $ 

DIAGONAL 0MEGAl/0MEGA2/*SQUARE*/2.0 ♦ 

MERGE MDD21rrMDD22r»VDf / MDD2 /1/0/2 $ 

MERGE BDD21* fBDD22»»VD* / BDD2 /1/0/2 ♦ 

MERGE KDD21»fKDD22**VD» / KDD2 /1/0/2 % 

SMPYAD MDD2fUDVFfOMEGA2*r» / TEMP2 / 3/-1 ♦ 

SMPYAD BDD2»UDVF»0MEGAI* » » / TEMP3 / 3 ♦ 

ADD TEMP2»TEMP3 / TEMP4 $ 

MPYAD KDD2*UDVF*TEMP4 / PPF2 / 0 ♦ 

MERGE PPFl*PPF2*r**VP / PPF / 1/0/2 $ 

EQUIO PPF*PNEF/MPCF1 ♦ 

COND LBLlSAfMPCFl ♦ 

DEC USETD/OECNEM/*P)tc/#NE*/*M* % 

PARTN PPF**UECNEM/PBNEF*PMF*r / 1 $ 

MPYAD GMD*PMF»PBNEF/ PNEF / 1 % 

LABEL LBL15A ♦ 

EQUiy PNEFrPFEF/SINGLE * 

COND LBLlSBfSINGLE $ 

VEC USETD/OECFES/*NE*/*FE*/*S# * 

PARTN PNEF* »VECFES/PFEF*PSF** / 1 $ 

LABEL LBL15B $ 

EQUiy PFEF*PDF/OMIT ♦ 

COND LBL15C.0MIT $ 

DEC IJSETD/ VECD0/»FE*/#D»/*0* $ 

PARTN PFEF* *OECDO/PBDF»POF» * / 1 $ 

MPYAD GOD*POF»PBDF/ PDF/ 1 $ 

LABEL LBL15C $ 

ALTER 101*101 $ USE FOL INSTEAD OF PPF 

ODR CASEXX * EQDYN * USETD * UHOF *FOL * XYCDB * /OUHUCl * /#FREORESP*/ 

♦MODAL*/S * N * N0S0RT2/S * N * NOH/S * N * NOP/FMODE ♦ 

ALTER 116*116 $ USE FOL INSTEAD OF PPF 

DDR2 USETD * UDUl F * PDF * K2DD * B2DD * MDD * FOL * LLL * DM/ 

UDU2F*UEUF*PAF/)|{FREQRESP#/N0UE/REACT/FRQSET $ 

ALTER 122*122 $ USE FOL INSTEAD OF PPF 

SDR2 CASEXX*CSTM*MPT*DIT*EQDYN*SILD* * *BGPDP*FOL*OPC*UPVC*EST* 

XYCDB*PPF/0PPC1 *00PC1 *0UPMC1 *0ESC1 *0EFC1 *PUGU/#FREQ#/ 
S*N*N0S0RT2 $ 

ALTER 128*128 $ 

SDR2 CASEXX*CSTM*MPT*DIT*EQDYN*SILD* * * *LAMA*QPH*PHIPH*EST»XYCDB*/ 
* IQPl * IPHIPl * * lEFl */*MMREIG»/S*N*N0S0RT2 ♦ 

ALTER 129*129 $ USE FOL INSTEAD OF PPF 

SDR2 CASEXX* * * *EQDYN*SILD* * * * FOL* * * * XYCDB* PPF/ 

OPPCA* * * * */#FREQ* $ 

VEC USETD / VD2 / *D*/#A*/#E* $ 
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VEC 

PARTN 

PARTN 

MPYAD 


USETlt / VA2 / *A*/*L*/»R* * 

PHIDHir »VIi2 / PHAVF»PHEVF»» / 1 $ 
PHAVFf»VA2 / PHL.VF»PHRVF»» /It 
DM»PHRVF>-PHLUF / PHLVFP / 0 / -1 t 


ORIGINAL PAGE ^ 
OF POOR QUALITY 


UMERGE USETD»PHLVFP» / PHAVFP / #A*/*L*/*R* t 
UMERGE USETDrPHAUFPfPHEyF / FHDVFP / *DJ|t/*AHt/#E* * 

SDRl USETD» fPHDUFP. » »G0D»6Ml>» » » » / PH1PHP» » / 1/*DYNAMICS» * 

SDR2 CASEXX » CSTM » MPT » D1 T » EQDYN» SILD » » » » I AMA » » PHI PHP » EST i. 

XYCriB»/»»»IESl»» /*MMREIG#/S»N»N0SGRI2 t 
ALTER 134fi34 t UsE FOL INSTEAD OF PPF 

DDRMM CASEXX f UHVF » FOL » IPHI P2 » IQP2 » IES2 , 1EF2 . X YCDD ,EST , MPT » D I T/ 

ZUP0C2»ZQPC2»ZESC2»ZEPC2» ♦ 

ALTER 138rl3d $ USE FOL INSTEAD OF PPF 

DDRMM CASEXX »UHi'F» FOL. IPHIPl » IQPl » lESl flEFl . .EST.MPT.DIT/ 
ZUPUCl .ZQPCl .ZESCl .ZEFCl. $ 


ENDALTER 
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$* ALTERS TO SOL i2 (NASTRAN APRIL 1982) FOR BASE ACCELERATION. 

«S 

ALTER 88 ♦ 

PARAM //*ADD*/KDEKA/N0UE/N0K2PP * 

PARAM //*ADn»/KDEK2/N0GENL/N0SIHP ♦ 

EQUIO <AA»KDD/KDEKA « 

ALTER 90*91 ♦ 

PARA.» //»AUD*/N0BGG/-l/0 ♦ 

PARA.i //*ADD»/N0K4GG/-l/0 $ 

GKAJl USETDf GM»GOf KAA» r MAAf f K2PPr H2PPf B2PP/ 

KDD » BDD » HDD » GHD » GOD » K2DD » H2DD » B2DD/*TRANRESP*/*DISP*/ 
*DIRECT*/C » Y f G=0 . 0/C » Y » W3=0 . 0/0 . 0/N0K2PP/N0H2PP/N0B2PP/ 
MPCF1/SINGLE/OMIT/NOUE/NOK4GG/NOBGG/KDEK2/V»Y»HODACC=-1 $ 
EQUIV B2DD»BDD/N0GPDT/H2DD»HDD/N0SIHP/K2DD»KDD/KDEK2 ♦ 

GKArt USETD » PHI A f H I f LAHA » D I T » » » » CASECC/ 

HXHH f BXHH » KXHHi* PHI DH/NOUE/C » Y » LH0DES=0/ 
C»Y»LFRE0=0.0/C»Y»HFREQ=-1.0/-1/-1/-1/ 
S>N»NONCUP/S»N»FHODE $ 

OEC USETD/VD/*D*/*COMP*/*R* 4 

EQUIV MAA»M1DD/N0UE ♦ 

COND LBL12rN0UE 4 

VEC iJSETD/VDl/»D»/#A»/*E* 4 

HERGE MAAj »»»VD 1,/ Ml DU /- 1/0/1 4 

LABEL LBL12 4 

PARTN M1DD»VD»/M1DD11» »M1DD12» 4 

MERGE M1DD11»»M1DD12» fODr/ M1DD2 /-l/O/l 4 

SMPYAD PHIDH!rMlDD2»PHIDH» » T/M1HH2/3/1/1/0/1 4 

DIAGONAL MXHH/MINV/)|cSQGARE#/-3 .0 4 

SMPYAD MIN0»BXHH»MlHH2f f »/BlHH/3 4 

ALTER *9f99 4 

TRLG CASEXXi.USrTDTDLT»SLT»BGPDT»SILrCSTM»TRLfDIT»GMDTGODf »ESTf 

MGG/PPT. ST»PDT»PD»»TOL / S»N»NOSET/NCOL 4 
ALTER 101 » 101 4 
PARTN PDf >0D/PD1»PD2» f/1 4 

PARTN MDDrVDf/MDDil»MDD21»MDD12»MDD22 4 

PARTN PDD»VDf/BDDll»BDD21fBDD12.BDD22 4 

PARTN KDD»0D»/KDDll»Ki:iD21 »KDD12»KDD22 4 

MERGE MDDllf f »MDD22»0Df / MDDl /-l/O/l 4 

MERGE fciDIill» »BDD12» f VDf/ BDDl /-l/O/l 4 

MERGE KDlill»»KD '^2»rOD»/ KDDI /-l/O/l 4 

MPYAD MDIH2»FD2»PD1/ PDDl /0/-1 4 

MPYAD MDD22 Rj2» / PDD2 /O 4 

MERGE PDD1»PDD2»»»»0D/ PDD /1/0/2 4 

MPYAD PK*DHfPDD»/PH/l/l/0/0 4 

SMPYAD r.,IDH)-MDDl»PHIDH»»»/MHH/3/l/l/0/l 4 

SMPYAD PHIDH»KDD1»PHIDH» » »/KHH/3/l/l/0/l 4 

EQUIV B1HM»BHH/N0B2PP 4 

CONI' LBL13A»N0B2?P 4 

SMPiAD PHIDHrBDDifPHIDHf » »BlHH/BHH/3/l/i/0/l 4 
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LABEL LBL13A $ 

PARAM //*ADB*/NONCUP/l/0 $ 

TRD CASEXX » TRL » NLET » DI T » KHH » BHHt MHH » PH/UHVT » PNLH/ 

♦MODAL #/NOUE/NONCUP/S » N p NCOL/C * Y » I START $ 

ALTER 128»128 $ 

SDR2 CASEXX»CSTMfMPTfDITfEQDYN»SILD,» »»LAMA»QPH)-PHIPH»EST» 

XYCDP»/»IQP1»IPHIP1»»IEF1,/*MMREIG* ♦ 

ALTER 129 ♦ 

VEC USETD / VD2 / ♦D»/*A*/*E* * 

VEC USETD / UA2 / *A*/*L*/*R» ♦ 

PARTN PHIDH»fVD2 / PHAUTf PHEVT» » /I * 

PARTN PHAVT.»VA2 / PHLVr»PHRUT» , / 1 ♦ 

MPYAD DM»PHRVT»PHLUT / PHLVTP / 0 / -1 ♦ 

UMERGE USETD»PHLOlP» / PHAVTP / *A»/»1 ♦/*R* $ 

UMEKGE USETD fPHAVTPiiPHEVT / PHDVTP / *D*/*A*/*E* ♦ 

SDRl USETD » »PHDUTP» » »GODfGMD» »» » / PHIPHPf. / 1/»DYNAMICS* ♦ 

SDR2 CASEXX»CSTMfMPTfDIT,EQDYN*SILD»»»»LAMAf »PHIPHP»EST» 

XYCDBf/»»»IESlf» /)!tMMREIG»/S»N»NOSORT2 $ 


ENDALTER 
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./ CHANGE NAME=FRi;i:?t< 

./ NUMBER SE0l--=70»NEUl = 71 ,INCR=1»INSERT=YES 
C 

C COMMON ERI.I2BC IS INITIALIZED BY ROUTINE FRRD2. 

C 

COMMON /FRD2BC/ IH 

./ NUMBER SEQ1=410»NEU1=411»INCR=1»INSERT=YES 
C 

C IF IHr^^Oi- USE COMPLEX DOUBLE PRECISION ARITHMETIC, 

C 

IF <IH ,EO. 0) ITY 4 

./ NUMBER SECU--=500rNEUl-510»INCR=l»INSERT=YES 
CALL MAKMCB ( MC » OUT » IROW r IFO » ITY ) 

*/ DELETE SE01==510»SE02=510 
./ CHANGE NAME=--FRD2C 

./ NUMBER SE01=120»NEW1=121»INCR=1»INSERT=YES 

C 

C COMMON FRD2BC IS INITIALIZED BY ROUTINE FRRD2, 

C 

COMMON /FRD2BC/ IH 

./ NUMBER SEQ.I- 24()!.NEWl=241»INCR=lfINSERT=YES 

C 

C IF IH=-0» USE COMPLEX DOUBLE PRECISION ARITHMETIC. 

C 

IF (IH .EO. 0) lOUT = 4 
./ NUMBER SEQl-37()rNEUl=371»INCR=l»INSERT=YES 

C 

C IF IH==0» THEN DO NOT USE INCORE CAPABILITIES BECAUSE FOR IH=0 

C COMPLEX DOUBLE PRECISION ARITHMETIC WILL BE REQUESTED 

C AND SUBROUTINE INCORE CALLED BY FRD2C IS WRITTEN ONLY 

C FOR COMPLEX SINGLE PRECISION MATRICES, 

C 

IF (IH ,EQ, 0) GO TO 102 
./ NUMBER BEQ1=1070»NEW1=1071»INCR=1»INSERT=YES 
C 

C IF IH^--()» USE COMPLEX DOUBLE PRECISION ARITHMETIC, 

C 

IF (IH ,EQ. 0) lOUT = 4 
,/ CHANGE NAME^FRRD2 

./ NUMBER SEQ1^27()»NEW1=271»INCR=1»INSERT=YES 
C 

C COMMON FRH2BC WILL BE USED BY ROUTINES FRD2B AND FRD2C. 

C 

COMMON /FRD2BC/ IH 

./ NUMBER SEQ:I”440vNEW1=441»INCR=1fINSERT=YES 
C 

C IF QHHL IS PURGED AND MACH NUMBER IS NEGATIVE THEN 

C SOLVE THE COUPLED EQUATION = (-W#»2,M + IW,B + K)U = P 

C USING COMPLEX DOUBLE PRECISION ARITHMETIC, 

C VARIABLE IH WILL BE USED TO CONTROL SOLUTION LOGIC IN 
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ORIGINAL PAGE 18 

ROUTINES FRRD2» ERD2B AND FRD2C. QUALITY 

IF <MCB<l).LE.O .AND. M.LT.0.0) NONCUP = 1 

NUMBER SE01=720»NEW1=721 » INCR=1 » INSERT=YES 

IF IH=0» THEN DO NOT USE INCORE CAPABILITIES BECAUSE FOR IH=0 
COMDEX DOUBLE PRECISION ARITHMETIC WILL BE REQUESTED 
AND SUBROUTINE INCORE CALLED BY FRD2C IS WRITTEN ONLY 
FOR COMPLEX SINGLE PRECISION MATRICES. 

IF (IH .EQ. 0) GO TO 20 
NUMBER SEQl=1210f NEW1=1211 » INCR=1 » INSERT=YES 

IF IH = 0* THEN CREATE NULL TRAILERS FOR SCR2<QHR> AND SCR3<QHI). 

THESE DATA BLOCKS ARE NORMALLY GENERATED BY FRD2A IF 
IH IS NOT EQUAL TO ZERO EACH TIME THRU THE LOOP ON 
NFREQ IN THIS ROUTINE. SINCE FRD2A IS NOT EXECUTED IF 
IH EQUALS ZERO THEN AFTER THE FIRST PASS THE TRAILERS 
FOR SCR2 AND SCR3 WOULD BE INCORRECT SINCE SCR2 AND 
SCR3 ARE ALSO USED BY FRD2C AS SCRATCH DATA SEfS. 

IF (IH .NE. 0) GO TO 38 
CALL MAKMCB(MCB»SCR2f0»0r0) 

CALL WRTTRL(MCB) 

MCB(l) = SCR3 
CALL WRTTRL(MCB) 

38 CONTINUE 

/ NUMBER SEQl=1480rNEWl=1481»INCR=l»INSERT=YES 

CREATE A PSEUDO FRL DATA BLOCK ON SCRl FROM DATA BLOCK FOL FOR 
INPUT TO ROUTINE FRD2F. (NO TRAILER IS NECESSARY). 

THE FREQUENCIES FROM FOL HAVE BEEN READ INTO IZ(1) DURING THE 
SET-UP AT THE BEGINING OF THIS ROUTINE. THESE FREQUENCIES MUST 
BE CONVERTED TO RADIAN FREQUENCIES FOR FRL (W = 2PI*F>. 

CALL G0PEN(SCR1 flZdBUFl) »1) 

DO 210 I = If NFREQ 
Z(I) = Z(I) * TWOPI 
210 CONTINUE 

CALL WRITE ( SCRl »Z»NFREQ»1) 

CALL CLOSE (SCR1»1) 

CALL FRD2F(MHH» BHHf KHHr SCRlf Ir NLOAD» NFREQ » PHF r UHVF) 

./ DELETE SEQ1=1490»SEQ2=1490 
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Young ' s Modulus “ 1E8 FL 

Poisson's Ratio ■ 0.4 

-4 2 

Material Density = 8E-4 FL T 



Base Acceleration Input for Figures 1 and 7 


Figure 2. Example 1 
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RFs 8 and 11 with Base 
Acceleration ALTERS 
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Figure 5. Load Required on Mass m. of Example 2 (Figure 3), Frequency Response 







